Effective temperatures from A5 to G5 spectral 
types, using Balmer line profiles 

C. van 't Veer-Menneret, C. Bentolila and D. Katz 

DASGAL, Observatoire de Paris-Meudon 
France 

Abstract. We show how previous works (Fuhrmann et al., van 't Veer-Menneret 
& Megessier) demonstrate the efficiency of the use of Balmer line profiles for 
effective temperature determination. In agreement with them, we insist on the 
physical interest of this method based on the behaviour of these lines with the 
variations of the parameters involved in the treatment of the convective trans- 
port. The comparison between Fuhrmann's results and ours, independently 
obtained, exhibits a quite good agreement. We show new results of effective 
temperature, gravity and metallicities for a few of our programme stars, rang- 
ing from solar to overabundant metallicities. 
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1. Introduction 

The first step, in a detailed abundance analysis, is to obtain an accurate effec- 
tive temperature. The quality of the other fundamental parameters, gravity and 
metallicity, will follow. Reliable fundamental parameters are also required for 
internal structure models used for computing diffusion processes or for astero- 
seismology. 

The wings of the first four members of the Balmer line series are very sen- 
sitive to effective temperature (hereafter Toff), and almost insensitive to grav- 
ity changes for T^^g less than 8500 K, with a limit around 5000 K, from where 
Stark broadening becomes inefficient. Recent works (Fuhrmann et al., 1993, 
1994; van 't Veer-Menneret and Megessier, 1996, hereafter VM) have shown 
that Balmer Line Profiles (hereafter BLP) are also sensitive to the temperature 
structure of the models used to interpret them. This means that BLPs depend 
on the treatment of the convection transport and on the metallicity entering the 
Opacity Distribution Functions (ODF). 

In section 2 we describe and comment on these previous works. In section 3 
we present our results, and compare them to Fuhrmann's ones. 

2. Previous Works 

In a series of papers, Fuhrmann et al. (1993, 1994) and Axer et al. (1994) have 
largely demonstrated the reliability of deriving T^s from the first 4 BLPs. Their 
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most important finding can be summarized in this way: in order to fit simul- 
taneously the 4 observed profiles with the ones computed using a single model, 
they had to lower the mixing-length free parameter a of the convection from 
1.5 down to 0.5, its usual values being in the range 1.25, 2.0. Their extensive 
analyses show that their finding is vahd for a large range of temperatures and 
metallicities, mainly for atmospheres cooler than 6500 K, with solar and defi- 
cient metal content. 

In fact the H„ profile was found to be insensitive to the choice of a, whereas 
the three other lines were found to be strongly affected, the line depth increasing 
with decreasing a. This behaviour is due to the fact that the profile is formed 
above the convective zone, and the other lines inside it, as it was demonstrated 
by Fuhrmann et al. (1993). We refer to Figures 3 and 4 of the same paper, 
for an extensive view of the effects on computed BLPs of model parameters 
Teff, gravity, metallicity and the mixing length to pressure scale-height ratio 
a = l/Hp. 

The main interesting consequence is that the BLPs are effective indicators 
of stellar atmosphere structure, as long as H„ and following lines are jointly 
interpreted, and constitute a direct evidence of the depth stratification of a 
stellar atmosphere. 

In a more recent paper Fuhrman et al. (1997) achieved these analyses by us- 
ing the Mg lb triplet for gravity dctcmination. Indeed these pressure-broadened 
lines are powerful gravity indicators, assuming that Tgff has been determined 
independently and that weak Mg i lines are available for Mg abundance deter- 
mination. 

In VM, with our own observational material and reduction procedures, we 
reached the same conclusions as Fuhrmann et al. (1993, 1994) about the in- 
fluence of model structure on BLPs and consequently on derived Teg, as for 
instance for the Sun and Procyon. We have confirmed the necessity to lower the 
value of a down to 0.5 for a simultaneous fitting of Hq, and H/3. We extend this 
result for hotter stellar atmospheres with overabundant metal content as well. 

3. Our Method and Results 
3.1. Observations and reduction 

The observational material was described in VM. We recall that we used the 
CCD receptor combined with the spectrograph AURELIE attached to the Coude 
Focus of the 152-cm refiector at the Observatoire de Haute-Provence. We ob- 
served in Hq, and ranges with a spectral resolution of 18 000, and in the 
613.5 nm and 552.0 nm ranges with a resolution of 30 000 for abundance anal- 
yses, the signal-to-noise ratio being around 400 for most of the exposures. 

We present in this paper some results for stars which were former targets for 
the unfortunate EVRIS space mission, and the analysis of which remains inter- 
esting to pursue. Some of them may be candidates for future asteroseismology 
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Table 1. The occurence of effects on BLPs of parameters involved in model compu- 
tation 
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programs, others for programs concerning 5 Scuti stars and chemically peculiar 
stars. 

We recall (see VM) that we used Kurucz' code ATLAS9 (1993) for comput- 
ing the models entering the calculation of profiles by the code BALMER9. As 
explained in VM, the study of convection treatment in ATLAS9 showed that 
the overshooting option used hypotheses which were not physically acceptable. 
Then we decided to remove this option in the computation of the models we 
needed. This difficulty was confirmed by Castelli et al. (1997). 

We mention as an information that we have carried out, under the UNIX 
system, an efficient automatisation of Kurucz' codes, allowing computations of 
grids of models, fluxes and Balmcr lines. 

3.2. Procedure 

Table summarizes the occurence of an effect on Balmer line broadening due 
to a change of the parameters involved in the model computation. We refer to 
Figure 4 in VM, which displays in what direction the variation of each param- 
eter acts, and how much. In Table ^ we have added 3 parameters entering the 
mixing length theory of convection: A/V is the surface over volume ratio of the 
convective element, y is a weighting factor on the optical depth in the convec- 
tion efficiency parameter averaged on optically thin and thick cases, and k is 
a factor on the convective flux taking into account the vertical motions of the 
elements. We refer to VM and Castelli et al. (1997) for more details, comments 
and references about them. They are free parameters, which act on the BLPs 
in the same way as a. In all cases presented here these parameters are chosen 
and flxed once and for all, only a was changed and the overshooting option was 
switched off. The metallicity and the gravity being chosen as starting values, 
we proceed as follows: (i) first of all Tcff is derived by fitting observed and com- 
puted Hq profiles, (ii) then we select the model which best fits both and 
by adjustment of a, (iii) the metallicity is derived by an iterative method 
through an abundance analysis which allows to determine the gravity and the 
microturbulence velocity. The results of such an iteration appear in the first 2 
rows of Table |[ where the 1st row gives the results before the iteration and the 
second one the results retained after the first and only iteration. 
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Table 2. Results for Procyon. The first 2 columns give the metallicity of the ODFs 
and the gravities used for the starting models, is the microturbulence velocity in 
kms~^. (*) refers to Fuhrmann et al. 1997 
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Table 3. Results for other programme stars. The second row gives the metaUicity of 
the ODF used, is the microturbulence velocity in kms^^. 
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3.3. Results 

Table || and || show respectively our results for Procyon, and for a few other 
stars from our programme. In all cases a = 0.5. Uncertainties can be summa- 
rized as follows, with [M/H] being the metallicity in the ODFs: A[A//iJ] = ±0.5 
leads to ATos = —200 K, and that leads to a decrease in the resulting [Fe/H] 
of -0.15 dex. A change in a of 0.75 corresponds to a change in T^s of about 
300 K, depending on Teg and metallicity. Figure 1 shows, in the case of Pro- 
cyon, the adequacy of the resulting model which represents both profiles. It is 
exactly as good for the other stars presented here. Moreover, we have to add 
that the ionization and excitation equilibria are well fulfilled, and that is the 
only disagreement with Fuhrman et al. (1997) which remains to explain. In the 
same way the discrepancy with the gravity deduced from the orbital parameters 
is not yet completely resolved. 

We recall that in VM we have shown the quite good agreement between 
the Tcff derived, for the 2 Am stars, from the 2 independent methods BLP and 
IRFM. The latter gave 7175 K for 63 Tau and 7025 K for r Uma, to be compared 
with the values from BLPs in Table ||. 
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Figure 1. Observed profiles (full lines) of Procyon fitted to the ones computed (dotted 
lines) with the 2nd model in Table 2 

4. Conclusion 

We conclude by stressing the remarkable internal consistency in the results of 
Fuhrmann and collaborators, likewise in our results, and the significant agree- 
ment between these two completely independant works. As a consequence, the 
value a = 0.5 for the mixing length parameter, can be considered as the recom- 
mended value to be introduced in the models used for T^g determination of all 
stars with Tgg between 8500 K and 5000 K, when Ha is not available. 
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